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ABSTRACT 


The lagged-averaged forecast (LAF) technique applied to tropical cyclone track 
prediction is a weighted sum of recent forecasts that were started from initial conditions 
at Various times lagging the start of the forecast period. The goal of this study is to re- 
duce the track forecast error at t +24 h. Two tests of the LAF approach are presented 
to demonstrate feasibilty. The first test uses the nine CLIPER forecasts initiated at 24, 
30, 36, 42, 48, 54, 60, 66 and 72 h prior to the common verifying time. In this test, the 
mean 24-h forecast error is reduced by 8 % relative to the 24-h CLIPER forecast above. 
In the second test, the “modified” LAF involves onlv the five CLIPER forecasts initiated 
at 24, 36, 48, 60 and 72 h prior to the verifying time. However, the 36-h through 72-h 
CLIPER forecasts are first modified using statistical regression equations that include 
predictors related to new track information since these forecasts were initiated. Signif- 
icant reductions in the track forecast error result from these statistical adjustments. The 
modified LAF applied to an independent sample results in a reduction from 189 km to 
124 km in the mean 24-h forecast error or a reduction of 34 %. This ts a significant 
improvement because the JTWC mean 24-h forecast error for the last four years is ap- 
proximately 210 km. The standard deviations are sigmificantly reduced from 118 km to 
69 km. Because the combination of the modified CLIPER forecasts in the LAF tech- 
nique results in a significant improvement in performance, it 1s recommended that this 


technique be applied operationally and also be tested with dynamical models. 
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I. INTRODUCTION 


Forecasters at the Joint Typhoon Warning Center (JTWC) on Guam have a difficult 
job in forecasting tropical cyclone movement. This job is made even harder by random 
and systematic errors in the model guidance used in preparing the forecast. Random 
errors can arise from the lack of adequate data to initialize the model whereas systematic 
error are caused by some intrinsic deficiency of the model (Peak and Elsberry 1982). 
Random errors can cause a ‘windshield wiper effect’ that rotates model forecast tracks 
alternately to the left and then to the nght of the actual path (Elsberry and Dobos 1989). 
The forecaster then is not certian if a forecast change in track orientation represents an 
actual turning or a Spurious oscillation. 

Several attempts have been made to reduce the systematic error in dynamical track 
prediction models. Renard ez al. (1973) used a post-processing technique to reduce svs- 
tematic errors in the Hurricane and Typhoon Track (HATRACK) forecasts. Elsberry 
and Frill (1980) used a series of regression equations to reduce the systematic error in 
the Fleet Numerical Oceanography Center (FNOC) Tropical Cyclone Model (TCM). 
Peak and Elsberry (1982) used a backward extrapolation technique that satisfactorily 
reduced systematic errors in the HATRACK, TCM and the Nested Two-way interactive 
Tropical Cyclone Model. 

The lagged-averaged forecast technique used with global numerical weather predic- 
tion models is one method of reducing the forecast error. Reduction of errors in the 
model guidance can aid the JTWC forecaster in his her quest to reduce track forecast 
errors. The removal of just the random error can save lives and millions of dollars by 
providing a more accurate and timely forecast of tropical storm movement. 

The main objective of this study is to use the lagged-averaged forecast technique to 
reduce random error in the western North Pacific CLImatology and PERsistance 
(CLIPER) model. A secondary objective is to provide the framework for future work 
in applving the lagged-averaged technique to a dynamical model, such as the One-way 
influence Tropical Cyclone Model (OTCM). 


Il. DESCRIPTION OF THE LAGGED-AVERAGED TECHNIQUE 


The lagged-averaged forecast approach was introduced by Lorenz (1977) who sta- 
tstically combined a series of previous forecasts to reduce errors. Hoffman and Kalnay 
(1983) proposed that a combination of forecasts verifying at the time of the present 
forecast might add valuable information to the present forecast. They called this ap- 
proach the lagged-averaged forecast (LAF) method. Each LAF ensemble member is an 
ordinary forecast started from initial conditions at a time lagging the start of the forecast 
period by a different amount. These forecasts are averaged at the common Verification 


time to obtain 
LAF = a;F; ale (1) 


In (1), @, is the weighting factor and F, is the corresponding LAF ensemble member. The 
LAF method is operationally feasible since the LAF ensemble members are produced 
during the normal operational cycle (Hoffman and Kalnay 1983). The LAF was tested 
as an alternate to Monte Carlo forecast (MCF) method. Jn the simplest form of the 
MCF method, an ensemble of initial states is selected randomly from a collection of 
possible initial states. Each ensemble member is then integrated in time and the ensemble 
of forecasts is used to calculate estimates of the desired statistics. This technique is very 
expensive because of the number of forecasts that have to be made. Tests by Hoffman 
and Kalnay showed that the forecast skill of the LAF is slightly superior to the MCF. 
Chen (1989) stated that the most difficult task in using the LAF is determining the 
best weighting factors for each of the past forecasts. For example, incorrect weighting 
factors mav degrade the accuracy of the LAF by weighting a forecast with a larger lag 
too much relative to a forecast with a smaller lag. The optimal weighting factors should 
result in a combination of the lagged-averaged forecasts that minimizes random errors. 
The lagged-averaged technique applied to tropical cyclone track prediction involves 
a series of recent forecasts that were started from initial conditions at various times lag- 
ging the start of the forecast period (Fig. 1). Consider a 24-h forecast from time t and 
all available forecasts during the past 48 hours that also are valid at time t + 24 h. 
During the past 48 h, eight other forecasts were generated that had a track position that 
also verified att + 24h. The lagged-averaged forecast for t + 24 h is the combination 


of the nine forecasts. The weights to be applied in (1) to each of the nine forecasts will 


be generated here through statistical regression methods. Notice also that these fore- 
casts can be validated at one or more past times from information available at time t. 
For example, the 72-h forecast made 48 h previously can be validated with known posi- 
tions at t - 36, t -24, t -12 and t = 0h (Fig. 1). These validations will be used in a 
modification of the tropical cvclone LAF. 
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well as the verifving times in relation to the t + 24h forecast. 


Ill. DATA BASE 


A. CLIPER FORECASTS 

The CLIPER technique is a prediction scheme based on a series of regression 
equations using predictors derived from persistence (past 12 h and 24 h storm move- 
ment), climatology (time of vear and storm location) and storm intensity (maximum 
sustained surface winds). The 24-h forecasts are based almost entirely on persistence, 
whereas climatology plays a bigger role at 48 h and almost completely dominates the 
technique at 72 h. 

CLIPER is used in this test of the feasibilitv of the LAF technique because it is a 
forecast technique that is available every six hours to the forecasters at JTWC. Conse- 
quently, nine CLIPER forecasts can be combined as in Fig. 1. Other techniques such 
as the dynamical models are available each 6 h, but they are based on new information 
only each 12 h. Only five such 12-h forecasts are then available for combination in the 
LAF technique. CLIPER may be considered as a no-skill type of forecast aid. If this 
lagged-averaged technique can improve CLIPER, then it is expected that the LAF could 
be applied to dynamically-based forecasts, such as those made by the the One-way in- 
fluence Tropical Cyclone Model (OTCM). 

Dr. T. Tsui and Mr. R. Miller of the Naval Environmental Prediction Research 
Facility (NEPRF) provided the files of warning positions used to generate the CLIPER 
forecasts. Mr. P.H. Dobos of the Naval Postgraduate School produced the forecasts 
using the western North Pacific CLIPER (or WPCLPR) technique developed by Xu and 
Neumann (1985). The CLIPER forecasts that are used in this study are generated using 
warning positions for the years 1984 - 87. A total of 109 storms occurred during this 
period (Table 1). 


Table 1. TROPICAL CYCLONES PER YEAR (ATCR 1984, 85, 86 AND 87). 
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The CLIPER forecasts are divided into dependent and independent data sets. The 
dependent set consists of all tropical cyclones that occurred during 1984 - 86. These 
storms accounted for 159§ CLIPER forecasts. Each forecast contains 12-, 24-, 36-, 48-, 
60- and 72-h forecast positions. In addition to the forecast positions, the data set in- 
cludes the current warning position and the warning positions for the previous 12 and 
24 h used to generate the CLIPER forecasts. The independent set contains the same in- 
formation for all tropical cyclones for 1987, which accounted for a total of 560 forecasts. 
The CLIPER forecasts in this study are based on warning positions in contrast to 
CLIPER forecasts used at JTWC, which use Combined Automatic Response to Query 
(CARQ) positions. The CLIPER generated with warning positions is not as accurate as 
the CARQ forecasts since the previous warning positions at 12 and 24 h are not updated 
with new information as is the CARQ. 


B. BEST TRACK DATA 

Dr. T. Tsui and Mr. R. Miller of NEPRF also provided the files of Best Track (BT) 
information.The BT positions are based on a post-storm analysis at JTWC that takes 
into account all of the available data, including aircraft fixes, satellite fixes, surface and 
aircraft radar fixes and even an occasional ship fix. The BT is a complete history of the 
storm track that represents only the large-scale motion because the smaller oscillations 
have been removed. 

The BT set includes 2308 positions in the dependent set and 1016 positions in the 
independent set. The reason for the larger number of positions in the BT set than in the 
CLIPER set is that a CLIPER forecast can be generated only after the storm has existed 
24h. 


C. MORECAST ERRORS 

The most widely used measure for verification of a tropical cyclone track is Forecast 
Error (FE), which is defined as the great circle distance between the forecast position and 
the best track position. In this feasibility study, FE is computed by taking the 24 h 
forecast position and subtracting it from the corresponding BT position (Fig. 2). The 
AX and AY are squared and added, and the square root of the sum is then the FE. The 
longitude (AX) position is corrected for latitude by multiplying by the cosine of the av- 
erage latitude between the BT position and the 24-h position. FE are compiled in terms 


of mean, median and standard deviation as measures of the improvement by the LAF 


method. The mean forecast error (MFE) is simply the total of all the FE divided by the 
number of forecasts. 

As discussed by Neumann and Pelissier (1981), FE is not an absolute measure of 
error, because the best track position contains uncertainties. As indicated above, the 
BT is a position determined during post-storm analysis at JTWC, and is subjectively 
smoothed to represent the overall large scale movement of the tropical cyclone (Sheets 
1986). Despite this limitation, FE will be used here as the measure of forecast skill. 
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Figure 2. Schematic diagram depicting forecast error. 


IV. PROCEDURES AND RESULTS 


A. LAGGED-AVERAGED FORECAST APPLICATION 

The first test of the lagged-averaged method uses a series of nine CLIPER forecasts 
verifving at 24h as in Fig. 1. For example, a 72-h forecast made at 0000 UTC on the 
17th will verify at 0000 UTC on the 20th. The 60-h forecast made at 1200 UTC on the 
17th will also verify at 0000 UTC on the 20th, along with the 48-h forecast made at 0000 
UTC on the 18th, the 36-h forecast made at 1200 UTC on the 18th and the 24-h forecast 
made at 0000 LTC on the 19th. 

In the first step, the CLIPER forecasts are arranged so that for each date-time group 
(DTG) the entry consists of the 24-, 36-, 48-, 60- and 72-h forecasts all valid at DTG 
+ 24h. To get the DTG for each forecast to the correct day, a calendar program ts used 
to add one day to the DTG. Linear interpolation between these forecasts is used to cal- 
culate 18-, 30-, 42-, S4- and 66-h forecasts. All longitudes computed by linear interpo- 
lation are corrected for latitude by using the cosine of the average latitude between the 
two interpolated points. Thus, the data base contained all nine forecasts valid at one 
time as in Fig. 1. One restriction with this method is that all the information in the first 
54h of any storm is never used because the statistical regression program uses only the 
DTG entries with all nine forecasts in the regression technique. For a DTG with all nine 
forecasts vahd at the same time, it is necessary to bypass the first 54 h. This restriction 
causes short-lived storms to be eliminated from the data set. 

Using a linear regression of this data set, equations simular to (1) are produced for 
predicting the lagged-averaged 24 h latitude (NLAT) and 24 h longitude (NLON). The 
equation for NLAT is provided all the nine forecast latitudes as predictors, and the 
predictand is the BT latitude for the verifying DTG. A similar procedure is used for 
NLON with nine forecast longitudes as predictors, and the predictand is the BT longi- 


tude. The resulting regression equations are 
NLAT = 0.75 LT24 + 0.33 LT30 — 0.03 LT42 — 0.04 LT54 — 0.019 L766 2) 
NLON = 1.26 LN30 — 0.25 LNS4 — 0.005 LN66. (3) 


In the above equations, LT24, LT30, LT42, LT54, LT66, LN30, LN54 and LN66 are the 
latitude (LT) and longitude (LN) of the forecasts initiated at the number of hours indi- 


cated. These variables are from the complete group of 18 forecast positions (latitudes 
and longitudes) from 24 to 72 h. The ’MINITAB RELEASE 5.1’ regression technique 
first checks the predictors for correlation with other predictors (Minitab 1985). This 
correlation is called “multicollinearity.” If the correlation is very high (an R-Squared 
value greater than 99.99%), the predictor is eliminated from the regression equation. 

The analysis of variance tables for each regression are given in Table 2 and 3. In 
the analysis of variance tables, the sequential sums of squares (SEQ SS) are shown. The 
first line gives SS(X2:X1), i.e., the reduction in the SS residual due to the fitting of the 
X2 term. The next line gives SS(X3:X1,X2), i.e.,the reduction in SS residual due to the 
fitting of the X3 term, given that X2 has already been added. The next line (if any) is 
SS(A4:X1,X2,43) and so on. 


Table 2, ANALYSIS OF VARIANCE FOR THE LATITUDE REGRESSION 
EQUATION, 
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Table 3. ANALYSIS OF VARIANCE FOR THE LONGITUDE REGRESSION 
EQUATION. 
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The 24-h forecast errors from applying (2) and (3) to the dependent data set are 
shown in Table 4. This forecast error is calculated by comparing the regression derived 
NLAT and NLON with the BT positions at the verifying time. Also shown are the 24-h 


forecast errors of the same data set without application of the regression equations. 


Table 4. 24-H FORECAST ERROR USING THE LAF TECHNIQUE WITH 
THE DEPENDENT DATA. 


Unmodified 24-h Forecast Error (km) Mean Median Standard 
Deviation 


—— AF 24-h Forecast Error (km) Mean Median Standard 
Deviation 






The regression equation applied to the dependent data set decreased the mean 24-h 
forecast error by 8 %o and the standard deviation of the 24-h forecast error by 6 %o. Since 
the decrease was small, the independent data set was not run and this method was 


modified, as described in the next section. 


B. MODIFIED LAGGED-AVERAGED APPROACH 
The second approach is a modification of the nine lagged-averaged forecast tech- 
nique in Fig. 1. Only the five CLIPER forecasts at 24, 36, 48, 60 and 72 h are used be- 


cause the interpolated positions at the intermediate times did not appear to add 


significant new information. The calendar program is used to add | day to the 24h 
DIG, 1.5 days to the 36 h DTG and so on until 3 days are added to the 72 h DIG. 
Using the calendar program in this method eliminated the undesirable feature of elimi- 
nating the first 54 h of every storm as in the first test. 
1. Modified CLIPER Forecasts 

In this second test, the 36-, 48-, 60- and 72-h CLIPER forecasts will be first 
modified to take into account additional information about the actual storm positions 
from that time until the time the 24-h CLIPER forecast is initiated (Fig. 1). For example, 
the 12-h position in the 36-h CLIPER forecast can now be verified from the warning 
position from which the 24-h CLIPER forecast is initiated. In principle, this present in- 
formation can be used to modify the 36-h CLIPER forecast prior to its use in the LAF. 
To derive the desired modifications, the difference between the 36-h CLIPER forecast 
position and the verifying BT position after 36 h is used as the predictand. In this data 
set, AY corresponds to the latitudinal displacement and AY corresponds to the longi- 
tudinal displacement. The AY difference is corrected for latitude by multiplication by the 


cosine of the average between the BI latitude and the 36-h forecast latitude (Tig. 3). 
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Figure 3. Schematic diagram of the longitudinal and the latitudinal error, AY and 
AY, of the 36-h CLIPER forecasts relative to the BT positions. 
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In addition to the new knowledge about actual storm positions early in each 
CLIPER forecast, other predictors such as past positions, past displacements, CLIPER 
displacements and deviations of the CLIPER forecast positions are used in the statistical 
modification (Fig. 4). These predictors are similar to those used by Peak and Elsberry 
(1982) in a post-processing technique for adjusting tropical cyclone tracks. 

The past-position predictors defined in Table 5 include the previous positions 
of the storm that are available at the time of the forecast, which establish a general track 
orientation and speed. For example, the 12 and 24h past positions are available for a 
36-h CLIPER forecast made 12 h previously. Thus, the t = 0h (current position), t - 
12h, t- 24h andt- 36h are available as possible predictors to modify the 36-h CLIPER 
forecast in Fig. 4. 

The past-displacement predictors in Table 5 are computed using the current 
storm position (time t) and subtracting the previous positions at 12. 24 and 36 h in the 
case of a 36-h CLIPER forecast made 12 h ago. The zonal and meridional differences 
for these prior times will generate six predictors. 

The CLIPER-displacement predictors in Table 4 are computed relative to the 
Warning position at the time of the forecast, which is t - 12 h in the case of the 36-h 
forecast in Fig. 4. This initial position is subtracted from the 12-h forecast position, the 
24-h forecast position and the 36-h forecast position. These zonal and meridional dif- 
ferences lead to six more predictors. 

As indicated above, the deviations of early portions of the CLIPER forecasts 
from the known warning positions prior to the current time are expected to be key pre- 
dictors. In the case of a 36-h forecast made 12 h ago, the 12-h CLIPER forecast posi- 
tion (valid at time t) is subtracted from the current position at time t. These zonal and 
meridional differences generate two additional predictors. 

A total of 22 predictors are available for the latitude and longitude regression 
equations for modifving the 36-h CLIPER forecast (Fig. 4 and Table 5). Data sets for 
modifving the CLIPER forecasts at 48, 60 and 72 h are generated in a similar fashion 
using the forecast data sets along with the BT set. If the DTG of each set matched, the 
BT position was subtracted from the corresponding forecast position to calculate the 
predictands. The predictors for the 48-, 60- and 72-h regression equations are calculated 
in a similar fashion. A total of 30 predictors are available for the 48 h (Fig. 5), 38 pre- 


dictors for the 60 h (Fig. 6) and 46 predictors for the 72 h (Fig. 7) regression equations. 
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More predictors are available for the longer forecast because two more predictors of 
each type occur. Then the forecast error was calculated using each computed AY and 
AY (see chapter on error statistics for description of forecast error). The data were then 


run through a statistical program to generate the mean forecast error and the standard 


deviation. These means and standard deviations wilt be used to itlustrate the improve- 
ments in the modified CLIPER forecasts at 36, 48, 60 and 72 h. 





Figure 4. Schematic diagram of the 36-h predictors relating the CLIPER forecasts, 


past positions and BT positions for modification of the 36-h CLIPER forecast. 





Figure 5. Schematic diagram of the 48-h predictors as in Figure 3, except for the 
modification of the 48-h CLIPER forecast. 
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Table 5. PREDICTORS AVAILABLE FOR THE REGRESSION EQUATIONS 
MODIFYING THE 36-H CLIPER FORECASTS. 


Predictor Type Predictor | Reference Postion Number of each type 
Name of predictor 
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Figure 6. Schematic diagram of the 60-h predictors as in Figure 3, except for the 
modification of the 60-h CLIPER forecast. 
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Figure 7. Schematic diagram of the 72-h predictors as in Figure 3, except for the 
modification of the 72-h CLIPER forecast. 


Regression equations for modifying the 36-, 48-, 60- and 72-h CLIPER forecasts 
are gencrated using a standard stepwise regression technique from “MINITAB VER- 
SION 5.1’. All of the predictors (latitudinal and longitudinal) are entered in the re- 
gressions for the latitudinal (A}) and for the longitudinal (AX) modifications. The 
predictors are entered (forward stepped) one at a time based on a T-ratio value greater 
than 2.0. As a predictor is added to the equation, the T-ratios for that predictor, as well 


as all of the predictors in the equation at that time, are computed. This allows for the 
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removal during each step of a predictor that drops below the 2.0 T-ratio criterion. In 
addition to the T-ratio criterion, the amount of variance (R-squared) accounted for by 
the addition of each predictor must have contributed at least a one percent increase. 

Once the regression equations are produced, each DTG of data run through the 
regression equations will produce a new AX and AY. Then in the case of the 36-h 
CLIPER forecast, the new AX and AY are first converted to degrees and added to their 
respective longitudinal and latitudinal 36-h CLIPER forecast positions. This procedure 
is applied to the 48-, 60- and 72-h CLIPER forecasts to arrive at the modified forecast 
position data sets. This was done with the dependent and independent data sets. 

Using the above criteria, the regression equations for modifying the latitudinal 


and longitudinal positions of the 36-h CLIPER forecast are 


MOD36LT = —1.18 PD36C + 4.54 PP36B + 0.4 CD36C + 4.88 PP36E 

— 509.2 (4) 
MOD36LN = — 1.61 PD36D + 0.46 CD36F — 0.29 CD36E + 0.4 PD36F 

+°7.1 PP36E — 121.03. (5) 


The predictors (see Table 5) in the regression equations are listed in the order in which 
they entered the equation. In all the regression equations modifying the CLIPER fore- 
casts, the past-displacement (PDtt) predictors are selected as the primary predictors. The 
predictor chosen least often was the deviation of early portions of the CLIPER from the 
known warning positions prior to the current ume. The Analysis of Variance (AOV) in 
Table 6 contains the T-ratio and the amount of the variance that is explained by the 
addition of each predictor to the regression equation. The percent of the variance ac- 
counted for by the latitude regression equation is larger than for the longitude regression 
equation even though the longitude equation contains more predictors. The latitude 
equation accounts for more variance than the longitude equation in all of the equation 


sets. This is probably because the latitude errors are smaller than the longitude errors. 


Table 6. ANALYSIS OF VARIANCE FOR 36-H REGRESSION EQUATIONS 
FOR MODIFIED LATITUDE AND LONGITUDE. 


36 h Predictors (LT) T-Ratio Percent of variance accounted for by 
each oe 
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Percent of variance accounted for by the regression 65.21 

equation 

36 h Predictors (LN) T-Ratio Percent of variance accounted for by 
each predictor 
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Percent of variance accounted for by the regression 54.84 


equation 
Table 7 shows the results of the adjustments to the 36-h CLIPER latitude 
(MOD36LT) in (4) and longitude (MOD36LN) in (5) being added to the 36-h CLIPER 


positions for the dependent and independent data sets. These modified positions verify 


















at the same time as the unmodified 24-h CLIPER positions. The forecast error statistics 
are then calculated comparing the unmodified with the modified 36-h CLIPER. The 
mean forecast error and the standard deviation of forecast error both displayed a sig- 
nificant reduction (Table 7). The dependent data set had a 32 % reduction and the in- 
dependent data set had a 23 % reduction in the mean forecast error, along with 31 % 


and 20 % reductions of the standard deviations of the forecast errors. 


Table 7. LONGITUDINAL (AX), LATITUDINAL (AY) AND TOTAL FORE- 
CAST ERRORS (KM) FOR 36-H UNMODIFIED AND MODIFIED. 
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che. 





The equations for the modification of the 48-h CLIPER forecasts are 


MODISLT = 0.125 PD4SC + 5.35 PPASB — 1.19 PDASE + 0.34 CD4I8G — 479.9 (6) 
MODASEN = 0.09 PDASD + 0.54 CD4S8SH — 0.42 CD4SG — 1.13 PDASF 
+ 0.47 PP36E — 54.3. (7) 


The AOV (Table 8) shows the T-ratio for the first predictor picked in both the latitude 
and longitude equations to be below the 2.0 criterion. These T-ratios were greatly re- 
duced by the addition of the last predictor in the table in each case, but the amount of 
explained variance was so large for each that they were kept in the regression equations. 
The overall percent of explained variance is larger than for the 36-h equations and again 


the explained variance was greater for the latitude than for the longitude. 


Table 8. ANALYSIS OF VARIANCE FOR 48-H REGRESSION EQUATIONS 
FOR MODIFIED LATITUDE AND LONGITUDE. 


each predictor 


Percent of variance accounted for by the regression GOA9 

equation 

48 h Predictors (LN) T-Ratio Percent of variance accounted for by 
each predictor 


Percent of variance accounted for by the regression 63.25 


equation 
Table 9 displays the results of applying MOD48LT in (6) and MOD4S8LN in (7) 
to the 48-h CLIPER independent and dependent data sets. After adding the adjustment 





























to the 48-h positions, the mean forecast error is reduced by 45 % for the dependent data 
and 36 % for the independent sample. The standard deviation of the forecast error is 
also reduced substantially with reductions of 38 % for the dependent data and 29 % for 
the independent sample. The errors for 48-h modified CLIPER, which verifies att + 24 
hin Fig. 1, are larger than the modified 36-h CLIPER, that also verifies at t + 24h. 


Table 9. LONGITUDINAL (AX), LATITUDINAL (AY) AND TOTAL FORE- 
CAST ERRORS (KM) FOR 48-H UNMODIFIED AND MODIFIED 
CLIPER. 
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Reduction in standard deviation of Forecast Error | 86 


The ‘regression equations for the modification of the 60-h CLIPER forecasts are 


MODOOLT = 0.47 PD6OOE + 5.11 PP6OD — 1.51 PD60G — 6.97 CD60K — 418.1 (8) 
MODOOLN = — 1.05 PD60F + 0.38 CD60L — 0.53 CD60K — 6.97 CD60B 
+ 16.0 PP6OM + 3.4 CD60D — 121.03. (9) 


In Table 10, the percent of explained variance again is greater for the latitude 
than for the longitude equation. In the latitude portion of the analysis of variance, the 
addition of the final predictor to the regression equation, CD60K, causes the percent of 
variance accounted for by the next to last predictor, PD60G, to fall below one. The 


predictor, PD60G, was kept in the equation since it still exhibited a large T-ratio. 


Table 10. ANALYSIS OF VARIANCE FOR 60-H REGRESSION EQUATIONS 
FOR MODIFIED LATITUDE AND LONGITUDE. 


60 h Predictors (LT) T-Ratio Percent of variance accounted for by 
each i 


Percent of variance accounted for by the regression 78.77 

equation 

60 h Predictors (LN) T-Ratio Percent of variance accounted for by 
each predictor 
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Percent of variance accounted for by the regression 65.88 
equation 


Table 11 shows the results of the adjustments to the 60-h CLIPER latitude 
(MOD6OLT) in (8) and longitude (MOD6OLN) in (9) for the dependent and independent 


samples. These modified 60-h forecasts , as did the unmodified counterparts, verify at the 



















same time as the unmodified 24-h CLIPER positions. The mean forecast error and the 
standard deviation of forecast error both displayed a significant reduction. The depend- 
ent data set had a 50 % reduction and the independent data set had a 34 % reduction 
in the mean forecast error along with 43 % and 26 % reductions of the standard devi- 


ation of the forecast error. 
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Table 11. LONGITUDINAL (AX), LATITUDINAL (AY) AND TOTAL FORE- 
CAST ERRORS (KM) FOR 60-H UNMODIFIED AND MODIFIED 
CLIPPER: 


60 h Dependent AX | AY | Forecast Modi- Modi- Modified 
Error fied AX | fied AY | Forecast Er- 
ror. 
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The regression equations for the modification of the 72-h CLIPER forecasts are 
MODIZLT = 0.82 PD72G + 5.1 PP72F — 1.82 PD72L + 0.55 CD72M — 397.8 (10) 
MODT2LN = 1.85 DEV72H — 6.24 CD72B — 0.8 CD72G + 0.83 PDT2N 

+ 182.36. (11) 


The regression equation for the 72-h longitude is the only one in which the modification 
of the CLIPER forecast included a Deviation from CLIPER as a predictor. It is even 
picked as the primarv predictor. Similar to the 60-h latitude regression equation, the 72-h 
latitude regression equation has a predictor with the amount of variance accounted for 
by the predictor, PD72L, of less than one. This predictor is again kept in the equation 
Since its T-ratio is greater than two. As with all the other regression equations, the lati- 


tude equation accounts for more variance than the longitude equation. 


Table 12. ANALYSIS OF VARIANCE FOR 72-H REGRESSION EQUATIONS 
FOR MODIFIED LATITUDE AND LONGITUDE. 


Percent of variance accounted for by 
each predictor 
TT Le cl 


Percent of variance accounted for by the regression 83.28 
equation 
72 h Predictors (LN) T-Ratio Percent of variance accounted for by 
each predictor 


Percent of variance accounted for by the regression 68.47 
equation 


Table 13 displays the results of applying MOD72LT in (10) and MOD72LN in 
(11) to the 72-h CLIPER independent and dependent samples. After applying the ad- 





























justment to the 72-h CLIPER positions, the mean forecast error is reduced by 52 % for 
the dependent data and 32 % for the independent sample. The standard deviation of the 
forecast crror is also reduced substantially with reductions of 49 % for the dependent 
data and 26 % for the independent sample. The errors for 72-h modified CLIPER are 
the largest of all the modified data sets, which all apply at t + 24 hin Fig. 1. 


Table 13. LONGITUDINAL (AN), LATITUDINAL (AY) AND TOTAL FORE- 
CAST ERRORS (KM) FOR 72-H UNMODIFIED AND MODIFIED 
CLIPE'E. 
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2. Modified LAF Forecast 
These modified CLIPER forecasts are combined with the unmodified 24 h 
CLIPER forecasts to produce the “modified” lagged-averaged forecast (LAF) equations. 
Stepwise regression is used with BT positions as the predictands and the modified 
CLIPER forecasts along with the unmodified 24-h forecasts as the predictors. The re- 


sulting equations are 
N24LT = 0.3 MOD36LT + 0.2 MODOOLT + 0.21 LT24 + 0.16 MOD72LT 

+ 0.15 MODA8LT + 0.202. (12) 
N24LN = 0.42 LN24 + 0.32 MOD36LN + 0.01 MODOOLN + 0.07 MOD7T2LN 


FOO DISLN.— 1-47, (13) 
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Equations (12) and (13) compare to equations (2) and (3) in that they are based on a 
lagged-averaged forecast technique, although (2) and (3) have nine available predictors 
and (12) and (13) draw on only five predictors. The latitude equation (11) uses 
MOD36LT as the number one predictor. This is something of a surprise since the LT24 
was expected to have a greater effect in the regression equation and it was picked only 
as the third predictor. In the longitude equation (13), LN24 was picked as the number 
one predictor as expected. 

In Table 14, the analvsis of variance for equations (12) and (13) is shown. The 
percent of variance accounted for by each of the additional predictors is less than one, 


but since the T-ratios are greater than two, the predictors are kept in the equations. 


Table 14. ANALYSIS OF VARIANCE FOR 24-H REGRESSION EQUATIONS 
FOR MODIFIED LAF LATITUDE AND LONGITUDE. 


each predictor 


Percent of variance accounted for by the regression 98.74 
equation 


to] inde 


= 


24h Predictors (LN) T-Ratio Percent of variance accounted for by 
each predictor 
156 
MODSGLN 
MODGOLN 
MOD72LN 
MODASLN 


Percent of variance accounted for by the regression 98.93 
equation 





Table 15 contains the results of applying (12) and (13) to the dependent and 
independent data. This resulted in a reduction of mean 24-h forecast error of 35 % for 
the dependent data and 34 %o for the independent sample. This is a significant reduction 


compared to the 8 % improvement from the first LAF, which was for the dependent 


24 


data. The reductions of standard deviation of forecast error are 34 % for the dependent 
data and 42 °% for the independent data, which may be compared to 6 % reduction for 
the onginal LAF. Thus, the modifications of the CLIPER forecasts used in the LAF 
technique result in a significant improvement in preformance. 


Table 15. 24H FORECAST ERROR USING THE MODIFIED LAF FORECAST 
TECHNIQUE. 


Unmodified 24 h Forecast Error Mean Median Standard 
Deviation 


independent ace | 


Modified LAF 24 h Forecast Error Mean Median Standard 
Deviation 


adependent Cid Cid 












Figures 8 and 9 display graphically the results of Table 15 for the dependent 
sample and the independent sample respectively. The solid line in each graph is the 
unmodified 24, 36, 48. 60 and 72-h mean CLIPER forecast errors. The dotted line re- 
presents the modified 36, 48, 60 and 72-h mean forecast errors. These two lines intersect 
since the 24-h mean forecast errors are not modified. A significant reduction in the mean 


24-h forecast errors by the application of the modified LAF is indicated in both samples. 


to 
Lo) 


=) 
2 
64 
{@) 
ne 
[S) 
{+ 
uw 
3 
fd 
oO 
ba 
4 
G 
3 


CLIPER FORECAST 
ERROR 


-0 0-000 e0 MODIFIED CLIFER 
FORECAST ERROR 


LAF 24-H FORE- 
CAST ERROR 


36 48 
FORECAST TIME (H) 





Figure 8. LAF 24-h error (star) for the dependent sample along with the mean 
forecast errors (km) for the unmodified CLIPER forecasts (solid line) and the mod- 
ihed 36-, 48-, 60- and 72-h CLIPER forecasts (dotted line). 
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Figure 9. LAF 24-h error as in Fig. 8, except for the independent sample. 
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Another indication of the improvement from the LAF technique is shown in 
Figure 10, which is a scatter diagram of the individual 24-h forecast errors for CLIPER 
and the “Modified” LAF for the independent sample. All the CLIPER forecast errors 
above and to the left of the diagonal line are improved by the application of the “Modi- 
fied’ LAF. The CLIPER forecasts errors below and to the right of the diagonal line are 
degraded. Although some small degradations are noted, the majority of the CLIPER 
forecast errors decreased after the application of the ‘Modified’ LAF. 
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Figure 10. Scatter diagram of the individual 24-h forecast errors for CLIPER versus 
the ‘Modified’ LAF for the independent sample. When two or more points have 
sinular values, the actual count is given. The lne indicates equal errors, and 


CLIPER forecast errors above the line are improved by the LAF technique. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The lagged-averaged forecast (LAF) technique applied to tropical cyclone track 
prediction involves a series of recent forecasts that were started from initial conditions 
at various times lagging the start of the forecast period. Consider a 24-h forecast from 
time t and all the available forecasts from the past 48 hours that are also valid at time t 
+ 24h. During the past 48 h, eight other forecasts were generated that had a track po- 
sition that also verified at the t + 24h. The lagged-averaged forecast fort + 24 his the 
combination of the nine forecasts. The weights applied to each of the nine forecasts are 
generated through statistical regression methods. The first test of the LAF technique 
with CLIPER forecasts resulted in an 8 % reduction in mean 24-h forecast error. 

To improve the LAF preformance, the CLIPER forecasts first are modified statis- 
tically to reduce the mean errors of the 36-, 48-, 60- and 72-h CLIPER forecasts valid 
at the initial time. The modified LAF combines these modified CLIPER forecasts along 
With the unmodified 24-h CLIPER forecasts in terms of a latitude and a longitude 
equation. These modified LAF equations applied to the dependent data set produce a 
reduction of the mean 24-h forecast error from 223 km to 146 km, or a reduction of 35 
“o. In addition, a reduction in the standard deviation from 152 km to 103 km indicates 
an improvement in consistency of the LAF forecasts versus the unmodified 24-h 
CLIPER forecasts. Although encouraging, these reductions are expected since the re- 
gression equations are derived using the dependent sample. The true test is the applica- 
tion to the independent sample, which results in a reduction from 189 km to 124 km in 
the mean 24-h forecast error and a reduction from 118 km to 69 km in the standard de- 
viation. This is a significant improvement because the JTWC mean 24-h forecast error 
for the last 4 years is approximately 213 km. The lower mean forecast error and stand- 
ard deviation in the independent sample from 1987 may be explained partly by the rela- 
tively easy storms, since the CLIPER performance was the best in a 10-year sample (R. 
Sheets, personal communication). Some of this improvement may be attributed to the 
six Supertyphoons during that year since CLIPER tends to perform well on 
Supertyphoons. However, the improvement is so large with the independent sample that 


the Modified LAF technique is concluded to have some ability to reduce forecast error. 


30 


B. RECOMMENDATIONS 

The Modified LAF technique in this study used CLIPER forecasts, which are a rel- 
atively no-skill tvpe of forecast. These CLIPER forecasts used warning positions instead 
of the CARQ positions used operationally at JTWC. Some improvement might be ex- 
pected from use of the operational CLIPER forecasts. However, this Modified LAF 
technique should be applied to the forecasts of a dvnamical model such as the One-way 
influence Tropical Cyclone Model (OTCM) to further improve the accuracy and con- 
sistency of the dynamical forecasts.. If longer track forecasts were available, the LAF 
technique could also be applied to 48-h and 72-h track forecasts that are even more im- 
portant for military evacuations and other preparedness actions in the western North 
Pacific area. 
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